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Abstract. Measurements of cosmic ray particles at energies above E 5 • 10 14 eV 
are performed by large area ground based air shower experiments. Only they provide 
the collection power required for obtaining sufficient statistics at the low flux levels 
involved. In this review we briefly outline the physics and astrophysics interests of 
such measurements and discuss in more detail various experimental techniques applied 
for reconstructing the energy and mass of the primary particles. These include surface 
arrays of particle detectors as well as observations of Cherenkov- and of fluorescence 
light. A large variety of air shower observables is then reconstructed from such 
data and used to infer the properties of the primary particles via comparisons to air 
shower simulations. Advantages, limitations, and systematic uncertainties of different 
approaches will be critically discussed. 



1. Introduction 

The cosmic ray (cr) energy spectrum extends from about 1 GeV to above 10 20 eV. Over 
this wide range of energies the intensity drops by more than 30 orders of magnitude. 
Despite the enormous dynamic range covered, the spectrum appears rather structureless 
and can be well approximated by broken power-laws dN/dE oc E^ 1 . Up to energies 
of some 10 14 eV the flux of particles is sufficiently high allowing measurements of their 
elemental distributions by high flying balloon or satellite experiments. Such studies 
have provided important implications for the origin and transport properties of CR's in 
the interstellar medium. Two prominent examples are ratios of secondary to primary 
elements, such as the B/C-ratio, which are used to extract the average amount of matter 
CR-particles have traversed from their sources to the solar system, or are radioactive 
isotopes, e.g. 10 Be or 26 Al, which carry information about the average 'age' of cosmic 
rays. With many new complex experiments taking data or starting up in the near future 
and with a possibly new generation of balloons, this is a vital field of research and has 
been subject to separate talks presented at this Symposium. 
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Above a few times 10 15 eV the flux drops to only one particle per m 2 per year. 
This excludes any type of 'direct observation' even in the near future, at least if high 
statistics is required. Ironically, one of the most prominent features of the CR energy 
spectrum is the steepening of the slope from 7 = 2.7 to 7 = 3.1 at an energy just above 
some 10 15 eV. This is known as the 'knee' of the spectrum. It was first deduced from 
observations of the shower size spectrum made by Kulikov and Khristianson et al in 
1956 but it still remains unclear as to what is the cause of this spectral steepening. 
At an energy above 10 18 eV the spectrum flattens again at what is called the 'ankle'. 
Data currently exist, though with very poor statistics, up to 3 • 10 20 eV and there seems 
to be no end to the energy spectrum |@. || . 

The origin and acceleration mechanism of these, so called, ultra- and extremely high 
energy cosmic rays have been subject to debate for several decades. Mainly for reasons 
of the required power the dominant acceleration sites are generally believed to be shocks 
associated with supernova remnants. Naturally, this leads to a power law spectrum as 
is observed experimentally. Detailed examination suggests that this process is limited 
to E/Z ^ 10 15 eV. Curiously, this coincides well with the knee at -Eknee — 4 ■ 10 15 eV, 
indicating that the feature may be related to the upper limit of acceleration. The 
underlying picture of particle acceleration in magnetic field irregularities in the vicinity 
of strong shocks suggests the maximum energies of different elements to scale with 
their rigidity R = pc/Ze. This naturally would lead to an overabundance of heavy 
elements above the knee, a prediction to be proven by experiments. A change in the 
CR propagation with decreasing galactic containment at higher energies has also been 
considered. This rising leakage results in a steepening of the CR energy spectrum 
and again would lead to a similar scaling with the rigidity of particles but would 
in addition predict anisotropics in the arrival direction of CR's with respect to the 
galactic plane. Besides such kind of 'conventional' source and propagation models 
@, H several other hypotheses have been discussed in the recent literature. These 
include the astrophysically motivated single source model of Erlykin and Wolfendale 
trying to explain possible structures around the knee, as well as several particle 
physics motivated scenarios which try to explain the knee due to different kinds of 
CR-interactions. For example, photodisintegration at the source |7j], interactions with 
gravitationally bound massive neutrinos ||, or sudden changes in the character of high- 
energy hadronic interactions during the development of extensive air shower (eas) [(J 
have been considered. 

To constrain the SN acceleration model from the other proposed mechanisms precise 
measurements of the primary energy spectrum and particularly of the mass composition 
as a function of energy are needed. Significant progress has been made on these problems 
in recent years, but the situation is far from being clear. 

The other target of great interest is the energy range around the Greisen-Zatsepin- 
Kuzmin (gzk) cut-off at E ~ 5 • 10 19 eV. Explanation of these particles requires the 
existence of extreme powerful sources within a distance of approximately 100 Mpc. Hot 
spots of radio galaxy lobes - if close enough - or topological defects from early epochs 
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of the universe would be potential candidates. This topic has been addressed at this 
Symposium in some detail by Ostrowski, Zavrtanik, and others. Therefore, I will be 
very brief on experimental aspects relevant to measurements in this energy range. 

2. Extensive Air Showers and Experimental Observables 

An air shower is a cascade of particles generated by the interaction of a single high energy 
primary cosmic ray particle with the atmosphere. The secondary particles produced in 
each collision - in case of a primary hadron mostly charged and neutral pions - may 
either decay or interact with another nucleus, thereby multiplying the particles within an 
eas. After reaching a maximum (X max ) in the number of secondary photons, electrons, 
muons, and hadrons, the shower attenuates as more and more particles fall below the 
threshold for further particle production. A disk of relativistic particles extended over 
an area with a diameter of some tens of meters at 10 14 eV to several kilometre at 10 20 
eV can then be observed at ground. This magnifying effect of the earth atmosphere 
allows to instrument only a very small portion of the EAS area and to still reconstruct 
the major properties of the primary particles. It is a lucky coincidence that at the 
energy where direct detection of CRs rays becomes impractical, the resulting air showers 
become big enough to be easily detectable at ground level. 

Due to the nature of the involved hadronic and electromagnetic interactions and the 
different decay properties of particles, an eas has three components, electromagnetic, 
muonic, and hadronic. On average, a 1 PeV primary proton will produce about 10 6 
particles at sea-level, 80% of which are photons, 18% electrons, 1.7% muons, and 
about 0.3% hadrons. Neutrinos will also be produced by weak decays of particles, but 
they remain unseen by ordinary EAS experiments. During their propagation through the 
atmosphere, relativistic charged particles will also produce Cherenkov light and, finally, 
excitations of the IP and liV-band of N 2 and N^, respectively, will give rise to emission 
of fluorescence light. 

Extracting the primary energy and mass from such measurements is not 
straightforward and a model must be adopted to relate the observed eas parameters 
to the properties of the primary particle. As shall be discussed below, the analysis is 
complicated by large fluctuations of eas observables and by the fact, that virtually all 
of the eas observables are sensitive to changes in the mass and energy of the primary 
particle so that a careful reconstruction of the energy spectrum requires knowledge of 
the chemical composition. 

Ignoring these complications for a moment, some basic characteristics of eas, as 
observed in Fig. ^ can be deduced from very simple considerations. On average, the first 
interaction of the primary particle with a nucleus in the atmosphere - mostly nitrogen 

- will occur within its characteristic interaction length, Xj oc 1/a. Due to the larger 
hadronic cross section of a Fe+N system as compared to p+N, the iron nucleus will 

- on average - interact higher in the atmosphere. Because of the strong absorption of 
the electromagnetic and hadronic shower components in the atmosphere, one expects 
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Figure 1. CORSIKA go) simulations of 
the longitudinal shower development using 
the QGJSET jll| hadronic interaction model. 
The two thick lines represent averages of the 
individual p and Fe showers. 



to detect fewer electrons and 
hadrons at ground in case of 
primary heavy nuclei. In a sim- 
plified picture, we furtheron 
may consider the primary iron 
nucleus as 56 nucleons each 
having l/5Q th of the primary 
energy. This has two im- 
portant consequences: i) the 
energy of produced particles 
within a single collision will be 
lower by about the same num- 
ber, so that more pions and 
kaons have a chance to decay 
into muons before reinteracting 
with another air nucleus, ii) the 
multiplicity of produced parti- 
cles within a single high energy 
collision roughly scales with 
the logarithm of the energy, 
so that 56 nucleons each with 
l/56 th of the energy, will pro- 
duce a higher number of sec- 
ondary particles, which again 
may decay into muons. As a 



result of the two effects one expects to observe more muons at ground in case of heavy 
primaries. Also, fluctuations in the number of particles will be smaller in case of heavy 
primary nuclei. Major experiments performing particle measurements at ground in the 
energy range of the knee include CASA-MIA [IT! 



EAS-TOP @, HEGRA fU|, KAS- 
CADE [[15], |i~6H , and Tibet-AS7 |HJ, only the latter two of which currently take data. 
Most of the data in the GZK energy range are based on the Akeno Giant- Air-Shower- 
Experiment (AG ASA) |Tj|. Their primary observables are lateral particle density dis- 
tributions at ground, p e ,fi,h( r ), and their integrals N e ^ h = f,^ 2nrp e ^ h (r)dr, yielding 
total (for extrapolations r\ — > and r2 — > oo) or truncated particle shower sizes. Trun- 
cated muon-sizes, A^ r , have been introduced by the KASCADE collaboration in order to 
avoid systematic uncertainties resulting from extrapolations of p M (r) to large distances 
not covered by the experiment. 

Besides measuring particles at ground, experiments also aim at observing traces 
of the longitudinal shower development. A basic parameter here is the position of the 
shower maximum, X max , which penetrates deeper into the atmosphere with increasing 
primary energy and decreasing primary mass. Experimentally, X max is often inferred 
from the lateral density distribution of Cherenkov photons, pch- Experiments following 
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Table 1. Compilation of experimental EAS observables and operating experiments. 
References are given in the text. 



Measurement of 


Observables 


Energy Ranges 


Operating Expts. 


charged particles 


shower size 

lateral densitv distr 

arrival time distr. 


£ 5 • 10 13 eV 


GRAPES, 

KASCADE 

Maket ANI, SPACE 


muons 


muon size 
lateral miion distr 
arrival time distr. 
/U-production height 


Z 5 • 10 13 eV 
£ 10 15 eV 


KASCADE, 
AGASA (> 10 18 eV) 


hadrons 


hadron size 
lateral hadron distr. 
energy distribution 
spatial correlations 


5 • 10 13 - 10 17 
(practical limit 
by detector 
area) 


KASCADE 


Cherenkov light 


Cherenkov size 
lateral Ch. distr. 
shower shape 


5 • 10 13 - 10 17 
(practical limit 
by area of Ch.-cone) 


BLANCA, HEGRA 
DICE, TUNKA 


Fluorescence light 


total light yield 
longitudinal profile 
time profile 


Z 10 17 eV 
(limited by 
signal/noise) 


HIRes 



Table 2. Compilation of EAS observables and their sensitivity to primary energy and 
mass. (W; energy estimates dependent also on primary mass) 



Energy determination Mass determination 



a-\gN^ + b-\gN e lgAyigiV; 

shape of electron lateral distribution; age-parameter 

mean muon production height 

muon arrival time distributions 

IgNh/lgNp, IgY^Eh/lgN^, spatial distr. of hadrons 

fluctuations of EAS parameters, e.g. lg(A^)/ ^g(A^) 
Pcti(120) ^ non-imaging counters: inner Cherenkov slope 

total Ch. light W Ch. telescopes: shape of shower image 

(the Ch. observables are often considered a measure of X n 
total fluorescence light position of shower maximum 



this approach include CASA-BLANCA [|TJ] and the HEGRA AIROBICC detectors J2D 
Also, air Cherenkov telescopes have been employed for sampling the shower maximum 
(HEGRA |fL4| , DICE [21]), and more recently, X max has been reconstructed from muon 
tracks at ground by means of triangulation (HEGRA p2j and KASCADE (talk presented 
by C. Biittner at this Symposium)). A very elegant approach is the detection of air 
fluorescence light at large distances from the shower axis giving true information about 
the longitudinal development of the EAS. The major drawback is the low light yield 
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limiting its application to energies E 10 17 eV. This technique pioneered by Fly's Eye 
is being employed by HIRes [E4j] and the Pierre Auger Observatory p5f 



A summary 

of the different experimental observables is given in table |l|. Deep underground muon 
detectors are omitted here because of limited space and because of their only indirect 
relation to eas observations. 



3. Techniques of Data Analysis and Selected Results 



3.1. Single parameter methods 

Until recently, most experiments have performed single parameter analyses using 
only one of the eas observables discussed above. The composition has then been 
inferred by comparing the mean value of the experimental observable to either different 
trial composition models or to mean values of simulated proton and iron primaries 
and interpolations to a hypothetical 'mean mass'. Obviously, a significant piece of 
information available from the full distribution of shower observables is ignored in 
such analyses. The problems here are manifold, most importantly: i) the extracted 
mean mass is ambiguous, i.e. the mean values of the experimental distribution can be 
described by a single 'mean mass' component, by a mix of pure protons and iron, or 
by any other combination. Clearly, a 'composition' in a strict sense is not obtained, 
ii) There is no clue of how well the simulations describe the experimental data. Only 
if the mean value of the experimental observable falls outside the window spanned by 

proton and iron simulations, doubts 
can be raised about the qual- 
ity of the simulations (or appa- 
ratus). Furthermore, publications 
discussing single parameter results 
often do another simplification by 
comparing measured and simulated 
observables not at the level of the 
detector, but instead comparing 
quantities which are inferred from 
experimental data by means of EAS 
simulations. A very prominent ex- 
ample of such kind of analyses are 
plots of X max as a function of en- 
ergy as shown in Fig. |^. In this and 
The mean depth of shower many other cases, X max is based 

nax & 

function of energy on measurements of the slope of 

deduced from CASA-BLANCA f§. The the Cherenkov lateral distribution 

lines represent the values for pure proton 
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Figure 2. 
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and pure iron samples as predicted by 
different hadronic interaction models. 



within about 120 m from the shower 
core. EAS simulations using differ- 
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ent primary particles are then performed to convert this experimental slope parameter to 
X max . Finally, this biased and model dependent X max -\ahxe is presented in such graphs 
and compared to X max -values obtained directly from CORSIKA simulations. Clearly, 
comparing experimental and simulated slope parameters directly instead of introducing 
secondary quantities appears much more appropriate and avoids sources of systematic 
biases. Determination of E prim in experiments using open Cherenkov counters is mostly 
based on the Cherenkov light intensity observed at about 120 m from the shower core 
where influences of the unknown primary mass are minimal. However, a careful analysis 
of this procedure exhibits a residual mass dependence on the order of ±10 % in the knee 
energy range [|27|- Also, uncertainties to the energy scale by at least the same amount 
are caused by hadronic interaction models f26| but, again, are not considered in plots 
of the type of Fig.^|. Finally, the figure also nicely demonstrates the aforementioned 
inabilities to judge the quality of the different hadronic interactions models or to infer 
a true composition. Despite this criticism, the qualitative feature of the data - a some- 
what lighter 'mean mass' towards the knee and a heavier one above - appears common 
to all interaction models. 

With increasing computer power, eas simulations have advanced in quality and 
quantity, i.e. both the complexity to which details of hadronic interactions and 
propagation effects in the atmosphere are taken into account has been improved as 
well as the number of simulated events available for comparison. Therefore, analyses 
of eas data could progress from comparing only mean values to full distributions of 
observables, the benefits of which have been discussed already. Two examples are 
presented in Fig. ^. The left hand side of the top panel presents the event-by-event 
distribution of the inner Cherenkov-slope (previously used as a measure for X max , see 
above) in an energy range 8 • 10 14 < E < 2 ■ 10 15 eV [2~B] and the right hand side 
the distribution of log / \ogN e obtained from KASCADE for zenith angles 18° - 25° 
and similar energies f28|. Both data sets are compared to CORSIKA simulations of 
different primary particles using the QGSJET model. Note, that the left hand tails of 
the distributions are well accounted for by proton simulations while the right hand tails 
are accounted for by iron simulations. The KASCADE distributions cover more than 
three orders of magnitude and thus provide some level of confidence in the results of the 
simulations. Both data sets also require at least 3 different primary masses to account 
for the full shape of the experimental distributions. However, it is also observed that 
the Cherenkov data alone do not provide sufficient discrimination power to distinguish 
the nitrogen from the iron distribution. The bottom panel shows the extracted mass 
fractions as a function of energy. Clearly, a heavier composition is observed at energies 
above the knee. Also, the relative abundance of iron is fairly similar in both data 
sets but the major difference is a reversed proton and Helium abundance resulting in 
a somewhat lighter composition from KASCADE data. To resolve this differences, 
several cross checks of the two data sets appear expedient, e.g. tests of the zenith angle 
dependence, still higher statistics in simulations to verify the tails of the distributions, 
etc. Such work is in progress now. 
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Figure 3. Event-by-event distributions of the inner Cherenkov slope (l.h.s.) [^6| and 
of the muon/electron ratio (r.h.s.) |^8). The lower panel shows the cumulative fractions 
of different primary particles as a function of their energy. 



3.2. Multi- component methods and non-parametric approaches 

Multi-component detector installations can measure several eas parameters simultane- 
ously on an event-by-event basis. Proper combinations can then be identified to provide 
an estimate of the primary energy and mass (see Table 1,2). The analysis of multivariate 
parameter distributions also needs to account for influences of fluctuations which may 
be different in different observables. Non-parametric Bayesian methods and neural net- 
work approaches are well suited for these purposes and they also specify the uncertainty 
of the results in a quantitative way [p9[| . Simpler approaches, like the KNN-method have 
also been employed fl30 |, but they suffer from various deficiencies like dependencies of 
the results on the density and thereby on the statistics of simulated data points. 

First comprehensive non-parametric approaches with promising results have been 
presented at this Symposium by M. Roth et al. In these techniques, each event 
is represented by an observation vector x = (Njf, N e , Nh . . .) T of n suitable EAS 
observables. This vector serves as input to a procedure based on Bayesian or neural 
network decision rules by which the observed event is assigned to a given class of 
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elemental groups, say p, 0, and Fe. The first step in such an analysis is a calculation 
of likelihood (probability density) distributions p(x\ui) by means of Monte Carlo 
simulations using different primary particles. Here, p(x\u>i) describes the probability 
to observe the multi-dimensional observation vector x in an event belonging to a given 
class Ui G {p, O, Fe}. The term 'non-parametric' indicates that the representations 
of the distributions (like probability density functions of Bayes classifiers or weights 
of neural networks) are no more specified by a-priori chosen functional forms, but are 
constructed through the analysis process by the given (simulated) data distributions 
themselves. Using the Bayes theorem one can then translate the likelihood p{x\uj,j) 
for finding an event x in a given class u>i to the actually required distribution p(ui\x) 
representing the probability of class Ui being associated to a measured event x. In this 
step, an assumption is made about an a priori knowledge of the relative abundance 
of each class, a subject being handled by Bayesian inference procedures. If there is 
no further knowledge, the prior probabilities should be assumed to be equal (Bayes' 
Postulate) . 

Analysing an experimentally measured event, a statistical decision on the primary 
particle type/energy is then to be made. The applied Bayesian approach of the 
statistical inference provides the Bayes optimal way of combining prior and experimental 
knowledge and the Bayes theorem specifies how such modification should be made 
In simplest case, the Bayes optimal decision rule is to classify x into class Ui, 
if p{uJi\x) > p{ujj\x) for all classes Uj ^ oji. In this way one is also able to specify 
the actual quality of the decision by classification and misclassification matrices of the 
results. Unfortunately, because of limited computer time, analyses making use of all 
these pieces of important information have not yet been performed for high statistics 
data. However, good agreement to faster neural network approaches was reported in 



Ref. [32 



It is important to realize that the composition extracted that way yields a somewhat 
heavier composition as compared to the results from Fig. [| This effect may possibly be 
caused by different treatments of correlations among the fluctuating eas observables. 
Furthermore, detailed investigations also show that within a given interaction model, 



different sets of EAS observables lead to different mass compositions | 32|| . Generally. 



inclusion of hadronic observables tends to shift the composition to a heavier one, a 
result already suggested by earlier studies of the KASCADE collaboration f33fl . These 
inconsistencies appear to be caused by deficiencies of the employed hadronic interaction 
models. Both kinds of effects are subject to further studies and, at present, any 
quantitative result on chemical composition may be subject to further changes. 

Results on chemical composition are often presented in terms of the mean 
logarithmic mass. However, this quantity has some severe drawbacks as any information 
about individual abundances is lost. Furthermore, in absence of an unique prescription 
on how to calculate the average from abundances of different mass groups, uncertainties 
in (In A) may be on the order of 0.2 and more. This has been source for some confusion 
in the past and can best be avoided by providing tables of reconstructed abundances of 
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A compilation of the all- 
particle energy spectrum re- 
constructed from various ex- 
periments is presented in fig- 
ure I]. The agreement appears 
reasonable and deviations are 
mostly explained by uncer- 
tainties in the energy scale by 
up to 25%, e.g. CASA MIA 
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Figure 4. Primary energy spectrum as obtained 
from different experiments fjifl . (See text for 
references.) 



data |35| were shifted upwards 
in energy by 20 % to yield a 
better agreement to the other 
data sets. This is likely to 
be explained by the interac- 
tion model SIBYLL 1.6 |§ 



employed by the authors of 



Ref. p5| but which has been proven to provide only a poor description of the exper- 



imental data ||33|| . The lines correspond to a simultaneous fit of the electron and muon 
size spectra of KASCADE, assuming the all-particle spectrum to be described by a sum 
of proton and iron primaries p3| . Interestingly, a knee is only reconstructed for the light 
component and no indication of a break is seen in the heavy component up to about 
10 17 eV. This important finding giving direct support to the picture of acceleration in 
magnetic fields (see above) will be target of future studies with improved experimental 
capabilities. For example, KASCADE and EAS-TOP have just started a common effort 
to install the EAS-TOP scintillators at the site of Forschungszentrum Karlsruhe pro- 
viding a 12 times larger acceptance as compared to the original KASCADE experiment 
and still taking advantage of the multi-detector capabilities. 



4. Summary and Outlook 

Much progress has been made in reconstructing the parameters of primary CR-particles 
from eas observables. The advent of multiparameter measurements and realistic EAS 
simulations now allow for much more detailed investigations taking account for eas 
fluctuations and detector effects. Different experiments agree fairly well on the all- 
particle energy spectrum but there is still some uncertainty in extracting the chemical 
composition. This may partly be explained by uncertainties of hadronic interaction 
models employed in eas simulations and there appears a demand for more tests of 
hadronic interaction models by means of accelerator and/or eas data. Other sources of 
systematic uncertainties appear to be caused by incomplete data analysis techniques and 
biases affecting mostly single parameter measurements. As discussions and presentations 
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at this Symposium have shown, there is well founded hope to solve these imperfections 
thereby answering the question about the sources of CRs and the origin of the knee 
already in the very near future. 
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